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A disorder parameter detecting dual superconductivty of the vacuum is used as a probe to char- 
acterize the confining properties of the phase diagram of two color QCD at finite temperature and 
density. We obtain evidence for the disappearing of dual superconductivity (deconfinement) in- 
duced by a finite density of baryonic matter, as well as for a coincidence of this phenomenon with 
the restoration of chiral symmetry both at zero and finite density. The saturation transition induced 
by Pauli blocking is studied as well, and a general warning is given about the possible effects that 
this unphysical transition could have on the study of the QCD phase diagram at strong values of 
the gauge coupling. 
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I. INTRODUCTION 

Color confinement emerges as an absolute property 
of strongly interacting matter from experimental facts, 
but is not yet fully understood starting from the QCD 
first principles. Lattice QCD simulations, however, have 
given some evidence about confinement and have even 
predicted the presence of a finite temperature transition 
to a deconfined state of matter, which is currently under 
investigation in Heavy Ion Collision Experiments. De- 
confinement, has seen by lattice simulations, appears to 
take place at the same temperature where other impor- 
tant physical phenomena happen, like for instance the 
restoration of chiral symmetry. 

In the present study we address the question regard- 
ing the fate of confining properties in presence of a finite 
density of baryonic matter: that is important in order 
to understand the structure of the QCD phase diagram 
both in the region of low densities and high temperatures, 
which is relevant for heavy ion experiments, and in the 
region of low temperatures and high densities, where the 
study of confining properties could help characterizing 
the nature of matter in compact astrophysical objects. In 
particular we are interested in understanding how decon- 
finement induced by a critical density of baryonic matter 
compares to what happens at zero density, and how it is 
related to other possible transitions in the QCD phase 
diagram. Various possibilities are open in principle, like 
for instance a non-coincidence of the deconfinement tran- 
sition with the restoration of chiral symmetry. 

The topic has been already studied in previous litera- 
ture. In particular results about deconfinement at high 
densities and low temperatures have been obtained in 
Ref. [l[, while the relation between deconfinement and 
the chiral transition has been investigated in Refs. 0, Q . 
However previous studies have been based on the study of 
the expectation value of the Polyakov loop as a probe for 
the deconfinement transition: while that can indeed be 
used as on order parameter for confinement in the pure 
gauge theory, being related to the spontaneous break- 
ing of the center symmetry, the same is not true in the 



theory with dynamical fermions, where center symme- 
try is explicitly broken. For that reason in the present 
study we look for different order parameters which are 
constructed in the framework of specific mechanisms of 
color confinement and which may be valid also in the full 
QCD theory. 

One appealing mechanism, among others, is that based 
on dual superconductivity of the QCD vacuum, which 
relates confinement to the breaking of an abelian dual 
symmetry induced by the condensation of magnetic 
monopoles [J, la, Ig]- The possibility to define disorder 
parameters in this scenario has been studied since a long 
time, one parameter has been developed by the Pisa 
group 0, HJ and consists in the expectation value of an 
operator which creates a magnetic monopole, (Ai) 1 : that 
has been shown to be a good parameter for color confine- 
ment both in pure Yang- Mills theories [l(| and in full 



QCD [Fy, 12 1; similar parameters have been developed 
elsewhere (IllllGJ]- 

Our plan is to use (M) as a parameter for characteriz- 
ing the confining properties of the various phases in the 
QCD phase diagram. Since numerical studies of QCD 
at finite density are notoriously difficult because of the 
sign problem, which makes usual importance sampling 
simulations unfeasible, in the present study we restrict 
ourselves to the theory with two colors, where that prob- 
lem is absent. In principle no significant differences are 
expected for the confining properties of the theory when 
going from N c = 2 to N c = 3, where N c is the number 
of colors: for that reason we believe that our study could 
be relevant also for real QCD. 

In Section |TT] we recall the general properties of lat- 
tice QCD at finite baryon density as well as the specific 
features of the two color model. In Section IIHI we re- 
view the definition of the disorder parameter (M) and 



1 We change the usual notation for the disorder parameter, (/i), 
in order to avoid confusion with the notation for the chemical 
potential. 
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develop our strategy to study its properties at finite den- 
sity. In Section IIVI we present our numerical results: in 
particular in Sect ion llV A I we present a study of the disor- 
der parameter at zero density; results at non-zero baryon 
density are given in Section IIVB1 a determination of the 
chiral transition line and its comparison with the decon- 
finement transition line are presented in Section [IV CI In 
Section TlV Dl the disorder parameter (Ai) and other ob- 
servables are used to study the nature of the unphysical 
transition to the saturation regime taking place at high 
values of the chemical potential and the possible influ- 
ence of saturation on the physical transition is discussed. 
In Section IIV El we discuss the relevance of our results 
for the low temperature region of the phase diagram. Fi- 
nally, in Section |Vl we draw our conclusions. A partial 
account of our results has been presented in Ref. [la |. 



expansion techniques [26|, l27j and non-relativistic expan- 
sions m n, m . 

The problem is absent in QCD with two colors, since 
the gauge group is real: indeed the fermion determi- 
nant, being expressible like any other gauge invariant 
observable in terms of traces over closed loops, is real 
as well, and numerical simulations are feasible. For 
this reason two-color QCD has been widely studied in 
the past as a laboratory for real QCD at finite den- 
sity [H, !, [H HI [H, H |M S HI . Despite some pecu- 
liar features of the model, like the fact that baryons and 
mesons are degenerate, one still expects to learn relevant 
information about specific questions, like for instance the 
fate of topology @ or confinement at finite density. 

III. THE DISORDER PARAMETER (M) 



II. QCD AT FINITE DENSITY AND THE 
TWO-COLOR MODEL 



We will consider a discretized lattice action for two- 
color QCD at finite chemical potential defined as follows: 



where Sg is the pure gauge Wilson action, 



□ 



(i) 



(2) 



the sum being over all plaquettes, while the fermion ma- 
trix is defined, in the case of standard staggered fermions, 
as 

1 3 

Here i and j refer to lattice sites, v is a unit vector on 
the lattice, ry 8ll/ are staggered phases and U are gauge 
link variables; a/i and am are respectively the chemi- 
cal potential and the quark mass in lattice units. The 
grand-canonical partition function can be written, after 
integrating out fermions, as: 



Z = J VUe~ SG[u] det M[U] 



(4) 



In ordinary QCD the fermion determinant is complex for 
generic values of the chemical potential, thus hindering 
the use of numerical Monte-Carlo simulations. Various 
possibilities have been explored to circumvent the prob- 
lcm, like for instance reweighting techniques fUEll, the 
use of an imaginary chemical potential either for ana- 
lytic continuation % [TJ, S3, El, [H H, Hlor for re- 
constructing the canonical partition function [251 ] , Taylor 



The magnetically charged operator M. (x , t) , whose ex- 
pectation value detects dual superconductivity, is defined 
in the continuum as the operator which creates a mag- 
netic monopole in x, t by shifting the quantum field by 
the classical vector potential of a monopole, and can 
be written (see Ref. Q for details) as 



M (x, t) = exp 



- / d 6 yE±{y,t)b±{y-x) 



(5) 



with the electric field E± (jj, t) being the momentum con- 
jugate to the quantum vector potential. Its expectation 
value, when discretized on the lattice, can be expressed 
as the ratio of two different partition functions, 



(M) = Z/Z, 



(6) 



where Z is the usual QCD partition function, while Z 
is obtained from Z by a change in the pure gauge action 
Sg — > Sg, consisting in the addition of the monopole field 
to the temporal plaquettes at a given timeslice where the 
monopole is created. 

Being expressed as the ratio of two different parti- 
tion functions, the numerical study of (A4) is a highly 
non-trivial task, since M. gets significant contributions 
only on those configurations having very small statis- 
tical weight. While numerical methods have been re- 
cently developed which permit a direct determination of 
(M) (38|, we shall not use them in the present study since 
they involve the combination of several different Monte 
Carlo simulations, a task which in presence of dynami- 
cal fermions could be unpractical. We will instead study, 
as usual, susceptibilities of the disorder parameter, from 
which the behaviour of (A 7 !) at the phase transition can 
be inferred. 

For instance, being interested in (M) as a function of 
/3, as for the /x = phase transition, one usually mea- 
sures 0, !, H 
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lnZ=(S) s ~(S)s (7) 
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where the subscript indicates the pure gauge action used 
for Monte Carlo sampling. The disorder parameter can 
be reconstructed from the susceptibility p, exploiting the 
fact that one has exactly (M) = 1 at j3 = 



{M)((3)=exp[ / p{P'W 



(8) 



In particular p ~ in the confined phase means {M) 
0, a sharp negative peak of p at the phase transition 
implies a sudden drop of (M) and p diverging in the 
thermodynamical limit in the deconfined phase means 
that (M) is exactly zero beyond the phase transition. 

Studying (Ai) as a function of (5 is what is usually done 
if interested in the fate of dual superconductivity as the 
temperature is increased. Indeed, in the Euclidean path 
integral formulation of QCD, the physical temperature is 
related to the inverse temporal extension, T = l/(L t a), 
where a is the lattice spacing which for an asymptotically 
free field theory is a decreasing function of the inverse 
gauge coupling /3. For that reason the inverse coupling /3 
is usually adopted in place of T when studying the QCD 
phase diagram, the latter being an increasing function of 
the former. 

At finite temperature and density we are interested 
in studying the behaviour of (M) in the two parame- 
ter space (/3,/t), where (i = ap is the chemical potential 
in lattice units. For that reason we introduce the new 
susceptibility 



p - 9 MM) - dlnZ dlnZ 
dp, dp dp 



= (N g )s-(N g ) s (9) 



where N q is the quark number operator, i.e. according 
to the definition of Z given in Eq. (HJ): 



(N g ) = 



( ~9A~ 



(10) 



(an additional factor 2 is actually needed for the case 
studied in the present paper, which deals with 8 staggered 
flavors, see Eq. ([12"]) '). The dependence of (M.) on the 
chemical potential p can then be reconstructed as follows: 

(M>GM) - (M)(/3,0)exp [ PD (p')dp' ) 



QCD phase diagram. Indeed, it is quite natural to as- 
sume that the gradient of the disorder parameter, 



(ii) 



be orthogonal, in the — p plane, to the critical line, 
whose slope is then equal to —po/p- In the following we 
shall directly check this property on our numerical data 
and also make use of it to obtain testable predictions. 



IV. NUMERICAL RESULTS 

In order to perform numerical simulations we have 
adopted the usual Hybrid Monte Carlo algorithm. The 
partition function in Eq. ([4]) can be rewritten, introduc- 
ing pseudo-fermionic fields as 



Z 



VUVQe 



-S g [U]-$*(M' t M)- 1 1> 



= J VUe~ s ° [u] (det M[U]f 



(12) 



In presence of a real chemical potential the usual even- 
odd factorization trick for reducing the number of flavors 
cannot be performed, so that Eq. (fl"2)) actually describes 
a theory with 8 (degenerate in the continuum limit) 
flavours. The standard exact 4> algorithm described in 
Ref. [11] has been used. 

We have performed simulations on lattices Ljj x L t with 
L t = 6 and different values of the spatial size ranging 
from L s = 8 to L s = 16. The bare quark mass has been 
fixed to am = 0.07. 

Simulations on the smallest lattice (L s = 8) have been 
performed on a PC farm, making use of a numerical code 
obtained by adapting the publicly available MILC code 
for two colors and for the inclusion of a finite chemical 
potential. Simulations on larger lattices have been per- 
formed instead on the INFN apeNEXT facility in Rome. 

The observables we look at are, apart from the sus- 
ceptibilities of the disorder operator introduced in Sec- 
tion lllli the average Polyakov loop, the average plaquette 
and the chiral condensate: 



(13) 



so that, if the starting point at p — is in the confined 
phase ((A4)(f3, 0) ^ 0), the behaviour expected for pd(p) 
in correspondence of a possible finite density deconfinc- 
ment transition will be the same shown by p across the 
finite temperature transition. 

Assuming the presence of a (pseudo)critical line in the 
T — p plane where the disorder parameter drops to zero 
and dual superconductivity disappears, the two suscep- 
tibilities p and po can be used not only to locate the 
position of the line, but also to compute its slope, thus 
providing a more comprehensive information about the 



&L t L\ ^ N e 



as well as their susceptibilities 



(14) 



(15) 



(16) 
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XL^Li ((L-(L)Y), 



XP=L 3 s L t ((P-(P)f) 



(17) 



(18) 



Notice that in the case of the chiral susceptibility we have 
explicitly considered only the disconnected contribution. 



A. The deconflning transition at zero chemical 
potential 

It is a well known fact that in ordinary full QCD at 
zero baryon density, chiral symmetry restoration takes 
place at the same critical temperature as deconfinement, 
with the latter identified with the disappearance of dual 
superconductivity [ll|, Q3 • We will check again this fact 
for the theory with two colors, since this will be an im- 
portant reference information for our following analysis 
at finite density. 
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FIG. 2: Chiral condensate and average Polykov loop as a func- 
tion of /3 measured on a 16 3 x 6 lattice at p = 0. The vertical 
line indicates the pseudo-critical coupling /3 C as determined 
from a best fit to the peak of the chiral susceptibility. The 
Polyakov loop has been multiplied by a factor 15 the better 
fit in the figure. 
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FIG. 1: Plaquette, Polyakov and chiral susceptibilities on the 
16 3 x 6 lattice at p = 0. The chiral and the plaquette sus- 
ceptibility have been respectively divided by a factor 10 and 
multiplied by a factor 4 in order to fit in the figure. Curves 
corresponding to best fits of the peak positions are superposed 
to the numerical data. 

We show in Fig. [1] the peaks of the three susceptibili- 
ties Xci XL and xp defined above, obtained on a 16 3 x 6 
lattice, together with curves corresponding to best fits 
to the location of their peaks. Our estimate for the lo- 
cation of the transition, obtained through a fit to the 
chiral susceptibility, is (3 C — 1.582(2), to be compared 
to those obtained by fitting the Polyakov loop suscep- 
tibility [Pl = 1-587(4)) and the plaquette susceptibility 
(f3p = 1.575(5)). A clear drop of the chiral condensate 
and a rise of the Polyakov loop are also observed at f3 c , 
as shown in Fig. O The dependence of /3 C on the spatial 



FIG. 3: p parameter as a function of f3 on various lattice 
sizes. The vertical line corresponds again to the location of 
the chiral transition. 

size is not significant, as can be appreciated from Ta- 
ble U where we report a summary of the pseudo-critical 
couplings (and chemical potentials) obtained from our 
simulations. 

Let us now consider the fate of dual superconductivity. 
In Fig. [3] we show the behaviour of the susceptibility p 
as a function of f3 for three different lattice sizes. A clear 
peak can be appreciated, which deepens when increasing 
the lattice size and whose location is clearly coincident 
with that of the chiral transition. Moreover it is also ap- 
parent from the figure that p is practically independent 
of the lattice size in the low coupling region, confirm- 
ing that (M) 7^ in the thermodynamical limit in that 
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FIG. 4: Behaviour of p in the strong coupling (/3 = 0.5) and 
in the weak coupling (/3 = 2.5) region as a function of L s for 
p, = and /i == 0.15. p stays constant and close to zero in the 
thermodynamical limit at strong coupling, while it diverges 
linearly with L s at weak coupling. 



phase, while p strongly depends on L s , and in particular 
is linear with it, as shown in Fig. 01 in the weak coupling 
region, showing that {M) is exactly equal to zero in the 
thermodynamical limit beyond the transition (magnetic 
charge superselection [40|). Therefore (3 C seems to sep- 
arate two phases characterized by a different realization 
of the U(l) magnetic symmetry. 

To better appreciate the coincidence of the chiral tran- 
sition with the disappearance of dual superconductivity, 
we have tried a finite size scaling (f.s.s.) analysis of the 
critical behaviour of (M) around the transition tempera- 
ture. We can assume for (M) the following f.s.s. ansatz: 



{M) = L s i cj>((p c -f3)L 1 J») 

from which it can be easily derived 

p = L\l v 4>{{j3 c - P)L\I») . 



(19) 



(20) 



We have checked this ansatz on our data, obtaining the 
best possible agreement for v ~ 0.63 and f3 c ~ 1.584: 
a reasonable scaling is obtained, with deviations ob- 
served on the smaller lattice. In particular we estimate 
j3 c = 1.584(2) in good agreement with the location of 
the chiral transition given above. The fitted critical in- 
dex v seems to indicate an Ising 3D critical behaviour, to 
be compared to that taking place in the quenched limit 
(3D Ising) and the renormalization group prediction for 
the critical behaviour in the chiral limit (first order [4l|). 
However a similar finite size scaling is not observed for 
the other susceptibilities and we believe that a definite 
answer about the universality class of the transition can- 
not be given in the present context, also due to the rel- 
atively small spatial volume used (our largest aspect ra- 
tio is slightly less than 3). A more careful investigation 



FIG. 5: Finite size scaling analysis of p around the transition. 
The best possible scaling is obtained for v ~ 0.63 and f3 c ~ 
1.584. 



should be performed and we consider the present analy- 
sis, as well as that presented later for the finite density 
case, as only aimed at a quantitative estimate of the crit- 
ical coupling where superconductivity disappears. 



B. The deconfining transition at non zero chemical 
potential 

The two susceptibilities p and pu permit to study (M) 
either as a function of temperature at a fixed value of the 
chemical potential /t, or as a function of p at fixed tem- 
perature. Both strategies can be used to investigate the 
fate of dual superconductivity in presence of a finite den- 
sity of baryonic matter: the first could be more effective 
at small chemical potentials, where the possible transi- 
tion line starting at p, = should be almost parallel to 
the p\ axis, the second could be more convenient at larger 
chemical potentials. Actually a proper combination of 
p and po could be used to study the behaviour of (M) 
along any given path in the (3 — p plane so that one could 
even choose an optimal combination corresponding to a 
relevant direction around a critical point: however we 
shall limit ourselves in the present context to the sim- 
pler cases of either fixed temperature or fixed chemical 
potential. The study at fixed temperature has a particu- 
lar interest, since it may show how the disappearance of 
confinement (dual superconductivity) can be induced by 
simply increasing the density of baryonic matter. 

We shall first consider the case of a fixed chemical po- 
tential, p = 0.15. In Fig. [6] we show the chiral suscepti- 
bility obtained on a 16 3 x 6 lattice and compared to the 
same quantity computed at p = 0. A clear shift of the 
pseudocritical coupling can be appreciated, in particu- 
lar we obtain f3 c (p = 0.15) = 1.568(2), showing that the 
(pseudo)critical temperature lowers as the chemical po- 
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FIG. 6: Chiral susceptibility on a 16 3 x 6 lattice as a function 
of /3 for various values of p. Dotted curves correspond to best 
fit to the peak values. 



tential is increased. Data for the susceptibility p on the 
same lattice are shown in Fig. and compared to those 
obtained at zero density: the peak of p shifts consistently 
by an amount comparable to that of the chiral suscep- 
tibility. Notice that in both cases the actual position of 
the p peak is at a (3 slightly larger than (3 C . That is ex- 
pected since p is a logarithmic derivative: assuming that 
(M)' = d(A4)/d(3 has a minimum at (3 C , it follows that 
dp/d(3 = (M)"/(M) - {{M)'/{M)) 2 is still negative at 
the same point. 

Data reported in Fig. [1] show that, also in the case 
fx = 0.15, p is independent of the lattice size and prac- 
tically equal to zero in the strong coupling region, while 
it diverges linearly with L s in the weak coupling region. 
Therefore we can conclude that, also in presence of a fi- 
nite density of baryonic matter, dual superconductivity 
disappears as the temperature is increased at the same 
point where chiral symmetry is restored. 

Next we turn to the behaviour of (M) as a function 
of p at fixed temperature (/?), determined by means of 
the susceptibility pn- We have considered only values 
of (3 below the (pseudo)critical coupling j3 c computed at 
jtt = 0, in particular (3 = 1.50 and (3 = 1.55: in this 
case we know that (M) at /t = 0, so that po may 
signal a possible disappearance of dual superconductivity 
induced by finite baryon density. Notice that the lowest 
value of /3, on the basis of a rough two-loop estimate 
of the /3-function, corresponds to a physical temperature 
T/T c ~ a(/3 = 1.582)/a(/3 = 1.5) - 0.4, where T c is the 
critical temperature at zero chemical potential. 

In Fig. [8] we show the chiral susceptibility determined 
on a 16 3 x 6 lattice at f) = 1.55 and (3 = 1.50. A 
best fit permits to locate the peak positions, hence 
the (pseudo)critical values of p corresponding to chiral 



restoration. We obtain p c ((3 = 1.50) = 0.340(10) 2 and 
p c {(3 = 1.55) = 0.215(10), as also reported in TableU 

In Fig. [3] we show instead the results obtained for pu 
as a function of p, at the same values of (3 and on various 
lattice sizes. It clearly appears that while po is inde- 
pendent of the lattice size and practically vanishing for 
small chemical potentials, it has a sharp negative peak 
in correspondence of the chiral transition which deep- 
ens as the spatial size is increased. In order to be more 
quantitative about the coincidence of chiral restoration 
and deconfinement, we have performed a f.s.s. analysis 
for the case (3 = 1.50, where three different spatial sizes 
were available (L s = 8, 12, 16), according to the ansatz 

(M) = i;^((Ac-A)^) (2i) 

hence 

PD = iy>((Ac - fi^vn . (22) 

A reasonable scaling is obtained for v ~ 0.55 and 
pL c ~ 0.31, in particular we estimate (3 C = 0.315(15), 
marginally compatible with the location of the chiral 
transition. 

We can therefore draw two important conclusions: 
dual superconductivity (confinement) disappears in pres- 
ence of a critical density of baryonic matter; moreover the 
critical line in the T—fj, plane corresponding to deconfine- 
ment coincides, at least within our present uncertainties, 
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FIG. 7: p on a 16 3 x 6 lattice as a function of /3 for various 
values of p. Vertical bands correspond to the pseudo-critical 
/3 fitted according to the chiral susceptibility. 



2 Notice that the (pseudo)critical chemical potential obtained at 
(3 = 1.50 is different from what obtained in Ref. |||: this differ- 
ence can be understood in terms of the different algorithm used, 
which in the case of Ref. was a non-exact molecular dynamics 
algorithm. 
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FIG. 8: Chiral susceptibility on a 16 3 x 6 lattice as a function 
of p for various values of /3. Dotted curves correspond to best 
fit to the peak values. 



FIG. 10: Finite size scaling analysis for pr>- A critical index 
v ~ 0.55 has been used, the best value for the critical chemical 
potential being f3 c ~ 0.315(15). 
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FIG. 9: as a function of fi and for two lattice sizes for 
various values of j3. Vertical bands correspond to the pseudo- 
critical chemical potential fitted according to the chiral sus- 
ceptibility. 



with the chiral transition line. These results concern that 
part of the phase diagram including temperatures down 
to T/T c ~ 0.4, where T c is the critical temperature at 
zero chemical potential: we shall discuss their relevance 
for the T ~ region of the phase diagram later in this 
paper. 



1.50) = 0.340(10), as obtained on our larger lattices (see 
Table IJ), we can perform a fit of the dependence /3 c (a0 
in the whole /3 — p. plane, which will then be used in the 
following. We are also interested in testing what stated 
in Section [ini i.e. that the ratio —po/p at the transition 
point can be used as an estimate of the slope of the criti- 
cal line: we give an example of a common plot of the two 
susceptibilites in Fig. [TTJ from which the ratio at (3 C , jj, c 
can be inferred. 
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FIG. 11: Comparison of p and pD as a function of fi at 
13 = 1.50 on a 16 3 x 6 lattice. 



C. The transition line 

Having obtained four different locations of the tran- 
sition line, in particular j3 c (p, = 0) = 1.582(2), (3 c (jl = 
0.15) = 1.568(2), p c (f3 = 1.55) = 0.215(10) and (ijfl = 



We have tried a quadratic fit /3 C (A) = A + Bfi 2 , obtain- 
ing A = 1.5828(16), B = -0.071(4) and x 2 /d.o.f. = 0.26. 
The good value of x 2 /d.o.f. shows that a quadratic de- 
pendence well describes the critical line down to T/T c ~ 
0.5; indeed a fit with a quartic term gives a coefficient 
for /t compatible with zero. Our estimates for the loca- 
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FIG. 12: Phase diagram in the /3 — p plane. The chiral line 
has been fitted to a quadratic dependence on fi. The slope 
of the critical line, has inferred from the disorder parameter 
for dual superconductivity, has been reported in the figure: 
a nice agreement (within one standard deviation) can be ap- 
preciated. 



tion of the (pseudo)critical points are reported in Fig.[T2l 
together with the fitted transition line. 

In correspondence of our direct locations of the tran- 
sition line we also show the estimates for the slope of 
the line obtained from the ratio — Pd/p' in particular we 
have drawn angles corresponding to one standard devia- 
tion from the average values. A good agreement can be 
appreciated, showing that —p^jp can indeed be taken 
as a good estimator of the slope of the line in the j3 — p 
plane. 

Finally, in Fig. 1131 we report again the chiral transi- 
tion line fitted above and compared to a quadratic fit in p 
for the critical line corresponding to the disappearance of 



CO. 



1.55 




1.45 



1.35^ 



FIG. 13: Comparison of the chiral pseudocritical line (contin- 
uous) and of that corresponding to the disappearance of dual 
superconductivity (dotted), as fitted from our data. 



dual superconductivity (deconfinement). The plot sup- 
ports our previous statement, i.e. that the chiral transi- 
tion coincides with deconfinement in the range of (3 values 
(temperatures) explored. 



D. A few remarks on saturation 

It is a well known fact that, even in absence of the sign 
problem, the study of lattice gauge theories in presence 
of a finite density of fermions cannot be pushed to ar- 
bitrarily high densities, i.e. to arbitrarily high values of 
the chemical potential. Indeed the number of available 
energy levels is limited by the presence of the UV cutoff, 
which places an upper limit to the possible values of the 
Fermi energy. Stated otherwise, we cannot place, because 
of the Pauli exclusion principle, more than one fermion 
with given quantum numbers per lattice site. Apart from 
the upper limit that this places on the densities reach- 
able on the lattice, a much worse problem comes from the 
fact that, as saturation sets in, the absence of available 
fermion levels quenches fermion dynamics, modifying the 
field theory at the ultraviolet scale. As a matter of fact, 
the theory becomes equivalent to a pure gauge theory in 
the large p limit. 

Saturation is therefore an unphysical lattice artifact 
which may in principle invalidate numerical results, one 
should therefore be extremely careful in locating its on- 
set. Indeed, while saturation effects are generically ex- 
pected to appear for p = ap of order 1, the exact value 
of p where they start to be important may depend on 
the dynamics of the theory. In the following we will 
briefly explore the transition to saturation in the two 
color model under consideration, arriving to some inter- 
esting conclusions which may sound as a general warning. 

We have explored saturation effects in some detail at 
(3 = 1.55. In Fig. [14] we show the behaviour of some 
observables as a function of p in a wide range going up 
to p = 1.6. For small values of the chemical potential 
above p c the fermion density rises roughly with a cubic 
dependence in p, as expected for a gas of free fermions, 
but then saturates to a value which in the figure is nor- 
malized to two fermions per site: the departure from the 
cubic behaviour starts at p ~ 0.6 — 0.8. Also the rise of 
the Polyakov loop suddenly stops at a similar value of 
p, followed by a drop; in the same region the plaquette 
suddenly drops towards its quenched value. Complete 
saturation is reached for p ~ 1.4 — 1.6. 

Much is learned by looking at the behaviour of the sus- 
ceptibilities of the disorder parameter in the same range, 
which is shown in Fig. 1151 the negative peak of po at 
p ~ 0.3, corresponding to the physical deconfinement 
transition, is followed by a positive unphysical peak at 
p ~ 0.7. That means that the disorder parameter (M), 
which at first drops to zero thus signalling deconfinement, 
then rises again as an effect of saturation: indeed the 
"saturation transition" leads to a SU(2) pure gauge the- 
ory, which at f3 = 1.55 and L t = 6 is deep in the confined 
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phase, implying (M) ^ 0. To verify that we have explic- 
itly reconstructed (M) (/t) / (M) (p, = 0) (see Ea.fTTJ) and 
reported it in Fig.[15] in the same figure we have reported 
the location of the saturation transition as obtained by 
a fit to the peak of the plaquette susceptibility. 

We should be satisfied, since the saturation transition 
at fX ~ 0.7 is well separated from the physical transition 
at fx ~ 0.3. However we notice that, defining a "satura- 
tion line" in the /3 — /t plane corresponding to the onset 
of saturation effects, we can predict, according to what 
stated in the previous paragraph, its slope from the ratio 
—po/p- We see from Fig. [T5]) that in correspondence of 
the positive saturation peak for pjj, the other suscepti- 
bility p has a negative peak, hence we expect a positive 
slope for the saturation line. That means that at lower 
values of (3 the onset of saturation could take place at 
lower values of ft: that combined with the fact that the 
physical critical f± c instead increases as (3 decreased, could 
lead to the unfortunate situation in which the two tran- 
sition, physical and unphysical, merge at lower values of 
(3, thus hindering, at least in the present case, the study 
of the strong coupling (low temperature) region of the 
phase diagram. 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 
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FIG. 15: p and pD at saturation. Notice the inversion of the 
peaks, indicating a transition line with positive slope. Ver- 
tical dotted bands refer to the locations of the physical and 
saturation transition. The thick continuous line refers to the 
disorder parameter (M) reconstructed by using the suscep- 
tibility po- after an intermediate region where the magnetic 
symmetry is restored, dual superconductivity sets in again in 
correspondence of saturation. 
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FIG. 14: Various observables computed at j3 = 1.55 on the 
8 3 x 6 lattice to show saturation. The Polyakov loop has been 
multiplied by a factor 8, the average plaquette by a factor 4. 

In order to further explore this possibility we have de- 
cided to make an estimate of the location of the satura- 
tion transition, through a fit to the plaquette suscepti- 
bilities (which are reported in Fig. I16j) . performing sim- 
ulations also at a different value of the gauge coupling, 
(3 = 1.675. Our estimate for the psudocritical saturation 
chemical potential /i£ c are reported in Table [J and are 
fi Sc ((3 = 1.55) = 0.68(3) and fi Sc (P = 1-675) = 0.79(3). 
In Fig. [T7] we report our estimate for the location of the 
saturation line together with a rough linear extrapolation 
suggesting that the saturation line could meet the phys- 
ical line, whose estimate given in previous paragraph is 
reported in the figure as well, for (3 ~ 1.4. Notice that the 



a n 

FIG. 16: Plaquette susceptibilities used to locate the satura- 
tion transition. 



linear extrapolation adopted is supported by the slope of 
the line obtained through the ratio —po/p, whose esti- 
mates are reported in the figure as well. 

We therefore give a general warning about the possible 
effects of saturation on the study of finite density QCD 
at low values of the gauge coupling. The situation may 
of course be quite different depending on the temporal 
extent L t of the lattice, on the number of flavors, of col- 
ors and on the lattice discretization (staggered or Wilson 
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FIG. 17: Saturation transition line and its relation with the 
physical transition line. 







fj, c 


8 3 x6 


1.584(2) 





12 a x 6 


1.587(2) 





16 a x 6 


1.582(2) 





16 a x 6 


1.568(2) 


0.15 


8 3 x6 


1.55 


0.222(10) 


16 a x 6 


1.55 


0.215(10) 


8 3 x6 


1.5 


0.325(10) 


12 a x 6 


1.5 


0.349(15) 


16 a x 6 


1.5 


0.342(10) 




/3 Sc 




8 3 x6 


1.55 


0.678(20) 


8 a x 6 


1.675 


0.793(10) 


16 a x 6 


1.675 


0.789(22) 



TABLE I: Collection of pseudocritical couplings as deter- 
mined from our numerical data. Physical critical couplings 
have been determined through the chiral susceptibility, while 
the unphysical saturation transitions have been located by 
means of the plaquette susceptibility. 

fermions) adopted. We plan to make a more extensive 
study of this problem in the future. 



E. Did we catch the physics of the low 
temperature region of the phase diagram? 

One of our starting questions was about the fate of con- 
finement at high densities and low temperatures, since 
that could help understanding the nature of compact as- 
trophysical objects. In Ref. [l| the hypothesis has been 
made, based on the analysis of the Polyakov loop, that 
at T ~ deconfinement could occur at a critical density 
following and well separated from the onset of a bosonic 
superfluid phase. It is natural to ask whether our present 



results can be of any relevance regarding this specific is- 
sue, i.e. how close we have got to the low temperature 
region of the QCD phase diagram. 

Since we have not included an explicit diquark source 
term in our model, we cannot obtain direct information 
about that observable; however we shall try to sketch a 
qualitative picture based on the distribution of the eigen- 
values of the fermionic matrix. At zero density that can 
be written as M — amid + D where D is antihermitian, 
hence it has purely imaginary eigenvalues, therefore the 
eigenvalues of M lie on a segment in the complex plane 
orthogonal to the real axis. 

As a real chemical potential is switched on, D ceases to 
be antihermitean and the eigenvalues get scattered in the 
whole complex plane: that is evident in the first inset of 
Fig. 1181 where we show the distribution of eigenvalues on 
a typical configuration obtained at (3 = 1.55 and a small 
chemical potential, £i — 0.10. The eigenvalues occupy a 
narrow vertical band and the finite density of eigenval- 
ues in correspondence of the real axis is strictly linked to 
the presence of chiral symmetry breaking (Banks-Casher 
relation (42J) The width of the distribution on the real 
axis grows as ft increases, roughly proportionally to /t 2 , 
till the distribution touches the imaginary axis: at this 
point the chiral condensate is expected to rotate into a 
diquark condensate (see for instance Ref. [43| for a re- 
view): as it is clear from the second inset in Fig. [TH1 
at (3 = 1.55 this happens roughly at /} ~ 0.3, a value 
which actually turns out to be almost independent of the 
gauge coupling in the range of (3 values explored in our 
simulations and is in agreement with the values found 
for diquark condensation in similar works using the same 
quark mass [3l|, |44|: we show as an example in Fig. [TH] 
the eigenvalue distribution projected onto the real axis 
for three values of jl at (3 = 1.45. 

However if we arc in the high temperature region, i.e. 
slightly below T c (fi — 0), chiral symmetry will be restored 
quite soon as /} is increased because of the transition to 
the Quark-Gluon Plasma. Therefore there will be actu- 
ally no chiral condensate to be rotated into a diquark 
condensate at the point where the distribution touches 
the imaginary axis. Indeed we see from the second inset 
in Fig. [TS] that the region around the real axis is quite 
depleted of eigenvalues for f3 — 1.55 at p, — 0.3. We can 
easily understand this in terms of the chiral line we have 
drawn in Fig. 1121 chiral symmetry gets restored already 
below jj, — 0.3 at j3 = 1.55. 

Following this line of reasoning, the region relevant 
for low temperature physics on our lattices with L t — 
6 should be that below (3 ~ 1.5, where our fitted 
(pseudo)critical line passes beyond /t ~ 0.3. In this region 
one could for instance observe, among other different pos- 
sibilities, two different transitions, the first corresponding 
to the onset of diquark condensation, the second roughly 
being the continuation of the line in Fig. Q~2J thus cor- 
responding to deconfinement: this is indeed the scenario 
suggested by Ref. [l[. 

We have therefore performed numerical simulations on 



11 



a 8 3 x 6 lattice at (3 = 1.45. In this case a finite density 
of eigenvalues around the real axis is still present at p ~ 
0.3, as can be better appreciated in Fig. [20j where we 
plot the distribution projected onto the imaginary axis 
at j3 = 1.45 and p — 0.3, compared to that obtained at 
higher temperatures. 

In Fig.[5T]we report the chiral susceptibility, compared 
to that measured on the same lattices at different gauge 
couplings. We notice that the peak is strongly reduced 
and its position is not much different from what obtained 
at (3 — 1.5 and in clear disagreement with what ex- 
pected from the continuation of the chiral line in Fig. [T2l 






FIG. 18: Distribution of eigenvalues on typical configurations 
of the 8 3 x 6 lattice obtained respectively at (/3 = 1.55, p = 
0.10), (/3 = 1.55, p = 0.30) and (/3 = 1.45, p = 0.30). 



The first peak could indeed correspond to the onset of 
a bosonic superfluid phase. Nothing seems to happen 
thereafter. 

In Fig.[2Uwe report instead data obtained for the sus- 
ceptibility pr> of the disorder parameter. In this case the 
negative peak has almost completely disappeared and a 
very small peak at p ~ 0.3 is followed by a region p > 0.4 
where po clearly changes its sign: on the basis of what 
we have discussed in Section IIVDI and comparing this 
behaviour with that observed at j3 = 1.55, a possible in- 
terpretation is that of an early onset of saturation effects 
in this case, preventing the observation of any further 
physical transition. We expected saturation effects to 
obscure the physical transition at (3 ~ 1.4, but we are 
not surprised that the situation may be worse. 

This conclusion is supported by looking at the be- 
haviour of the Polyakov loop (see Fig. |53J), in this case 
saturation effects are signalled by an inversion in the 
growth of (L) as a function of p. 
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FIG. 19: Distribution of eigenvalues projected onto the real 
axis for various p on a 8 3 x 6 lattice at j3 = 1.45. 
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FIG. 20: Distribution of eigenvalues projected onto the imag- 
inary axis for various values of (3 on a 8 X 6 lattice at p = 0.3. 
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We conclude therefore that we are not be able to clarify 
the onset of deconfinement at T ~ 0, at least on the 
present lattice size. We could of course further decrease 
the temperature without decreasing (3 by going to larger 
values of L s . However that would imply a numerical 
effort which is not affordable with our present algorithmic 
and computational resources. 



V. CONCLUSIONS 

We have investigated the phase diagram of two-color 
QCD at finite temperature and density by means of a 
disorder parameter for color confinement detecting dual 
superconductivity of the QCD vacuum. 

We have obtained evidence for deconfinement induced 




FIG. 21: Chiral susceptibilities on the 8x6 lattice at various 
values of 13. 




by a finite density of baryonic matter. Moreover the tran- 
sition line corresponding to the disappearance of dual 
superconductivity (deconfinement) appears to coincide, 
in the range of temperature explored (0.4 T c < T < T c , 
where T c is the critical temperature at zero density) , with 
that corresponding to chiral symmetry restoration, as it 
happens in the zero density case. We have also shown 
that the susceptibilities of the disorder parameter can be 
used in order to compute the slope of the critical line in 
the — jl plane, obtaining consistent results. 

We have investigated in some detail the unphysical 
transition corresponding to the onset of saturation and 
shown that it moves at lower values of jl as (3 is decreased 
with a possible intersection with the physical transition 
line, thus giving a general warning about the possible ef- 
fects of saturation on the study of finite density QCD at 
strong values of the gauge coupling. This phenomenon of 
course may be quite different depending on the fermion 
discretization, on the number of flavors and on other pa- 
rameters of the system (L t , quark masses) ; for this reason 
we plan to make a more systematic study in the future. 

We have also verified that in our case saturation ac- 
tually prevents us from obtaining results relevant for the 
T ~ region of the phase diagram. For this reason we 
plan to extend our study in the future by adopting differ- 
ent lattice sizes and/or fermion discretizations. Results 
relevant for the high temperature region of real QCD 
could also be obtained within the imaginary chemical po- 
tential approach. 
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